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Abstract: We continue an effort to obtain information on the QED per- 
turbation series at high loop orders, and particularly on the issue of large 
cancellations inside gauge invariant classes of graphs, using the example of 
the I - loop N - photon amplitudes in the limit of large photons numbers 
and low photon energies. As was previously shown, high-order information 
on these amplitudes can be obtained from a nonperturbative formula, due 
to Affleck et al., for the imaginary part of the QED effective lagrangian 
in a constant field. The procedure uses Borel analysis and leads, under 
some plausible assumptions, to a number of nontrivial predictions already 
at the three-loop level. Their direct verification would require a calculation 
of this 'Euler-Heisenberg lagrangian' at three-loops, which seems presently 
out of reach. Motivated by previous work by Dunne and Krasnansky on 
Euler-Heisenberg lagrangians in various dimensions, in the present work we 
initiate a new line of attack on this problem by deriving and proving the 
analogous predictions in the simpler setting of 1+1 dimensional QED. In the 
first part of this series, we obtain a generalization of the formula of Affleck 
et al. to this case, and show that, for both Scalar and Spinor QED, it cor- 
rectly predicts the leading asymptotic behaviour of the weak field expansion 
coefficients of the two loop Euler-Heisenberg lagrangians. 



1 Introduction 



QED is the oldest and prototypical quantum field theory, and much effort has 
gone into calculating low order terms in its perturbation series. However, 
little is still known about its high-order behaviour, in fact less than for 
some other field theories [TJ [3J HI [5] . This is due, on one hand, to the 
absence of (spacetime) instantons in the abelian case, on the other hand 
to large-scale cancellations between Feynman diagrams which are familiar 
to all practicioners in the field, but whose origin and effect on the large 
order behaviour of perturbation theory is still not well-understood. This 
type of cancellations has recently attracted wider attention since it is now 
increasingly believed that their study may throw light on the origin of similar 
cancellations observed in gravity and supergravity theories (see, e.g., [U [7] 
and refs. therein). 

In their pioneering calculation of the g — 2 factor of the electron to sixth 
order in 1974, Cvitanovic and Kinoshita [8] found a coefficient which was 
much smaller numerically than had been expected by a naive estimate based 
on the number of Feynman diagrams involved. A detailed analysis revealed 
extensive cancellations inside gauge invariant classes of diagrams. This led 
Cvitanovic |9] to conjecture that, at least in the quenched approximation (i.e. 
excluding diagrams involving virtual fermions) these cancellations would be 
important enough numerically to render this series convergent for the g — 2 
factor (see [10] for an amusing account of the genesis of this conjecture). 
Although nowadays there exist many good arguments against convergence 
of the QED perturbation series, all of them are based on the presence of 
an unlimited number of virtual fermions, so that Cvitanovic's conjecture 
for the quenched contribution is still open today. Moreover, should it hold 
true for the case of the g — 2 factor, it is natural to assume that it extends 
to arbitrary QED amplitudes, such as the quenched photon S-matrix (for 
photon amplitudes, we call 'quenched' the contributions involving only one 
scalar /fermion loop). Beyond the quenched approximation, it would suggest 
that the QED perturbation series should be rearranged as a series in the 
number of fermion loops, rather than in the coupling constant. 

The present work continues an effort [11] \12\ [131 fT4] to study the multi- 
loop behaviour of the QED N photon amplitudes using the QED effective 
lagrangian, and in particular to prove or disprove Cvitanovic's conjecture 
for these amplitudes. Let us start with recalling the representation obtained 
by Heisenberg and Euler [15] for the one-loop QED effective Lagrangian in 
a constant field, 
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Here T is the proper-time of the loop particle and a, b are defined by a 2 —b 2 = 
B 2 — E 2 , ab = E • B. The Euler-Heisenberg Lagrangian ('EHL' in the 
following) contains the information on the one-loop ./V photon amplitudes 
in the limit where all photon energies are small compared to the electron 
mass (see, e.g., [I6j [17]). After expanding the EHL in powers of the field 
invariants, the 'weak field expansion', it is straightforward to obtain the 
explicit form of the N photon amplitudes in this 'EH limit' from the terms 
in this expansion involving N powers of the field. Using the spinor helicity 
formalism, the result of this procedure can be expressed quite concisely [13]: 
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Here the superscripts ± refer to circular polarizations, and the are 
Bernoulli numbers. The invariants Xk are written, in standard spinor he- 
licity notation (see, e.g., [18j), 



f-V r i 
x + = 1^1|[12] 2 [34] 2 ---[(K-1)K] 2 + all permutations | 



2T 

X N -k = y -^&{{(K + l)(K + 2)) 2 {(K + ^){K + A)) 2 ...{{N-l)N) 2 
2~^~ 



+ all perm. | 

(1.3) 

The formula analogous to (jl.ip for Scalar QED is 
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This formula was obtained by Weisskopf [19|, but for simplicity will be called 
"scalar EHL" in the following. 

Except for the purely magnetic case, the effective Lagrangians (|l.ip . (|1.4p 
have an imaginary part. For the purely electric case, Schwinger [20] found 
the following representation of the imaginary parts in terms of infinite series 
of 'Schwinger exponentials', 
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These formulas imply that any constant electric field will lead to a certain 
probability for electron-positron pair creation from vacuum. The inverse 
exponential dependence on the field suggests to think of this as a tunneling 
process in which virtual pairs draw enough energy from the field to turn 
real, as had been proposed by Sauter as early as 1931 [21]. 

At the one-loop level, obtaining the imaginary part of the EHL from its 
real part is a simple application of Cauchy's theorem. For our multiloop 
purposes, however, it will be important that the imaginary part can, using 
Borel summation, also be obtained from the asymptotic behaviour of the 
coefficients of the weak-field expansion. Referring to [5j [22] for the details 
of this procedure, let us just cite the following main result: Assume that a 
function f(g) has an asymptotic series expansion 



(1.7) 



n=0 



where the expansion coefficients c n have the leading-order large n behaviour 



c n ~ p n T(pn + v) (1.8) 

with some real constants p > 0, p > and v. Then the series (II. 7p has 
convergence radius zero, and is not even Borel-summable. Nevertheless, one 
can, by a dispersion relation applied to the formal Borel integral, show that 
the leading imaginary contribution for small g is given by 
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These formulas can be applied to the weak field expansions of the purely 
electric EHL's (fTTTjl . (TQj) . setting g = (eE/m 2 ) 2 . The expansion coefficients 
c n involve essentially only the Bernoulli numbers [17} I22j . and using the 
asymptotic properties of those numbers it is easy to see that they obey 
(jl.8|) with p = l/V 2 , [i = 2, and v = —2, independently of spin. Eq. 
(jl.9p then reproduces the k = 1 terms in (|1.5p and (jl.6p . One can then 
go on and iterate this procedure to successively construct all the Schwinger 
exponentials [22J. However, in the following we will generally be concerned 
only with the leading term. 

Considerable work has gone into generalizing the EHL's to the two loop 
level, that is, taking a single photon exchange in the loop into account. For a 
general constant field, various equivalent integral representations have been 
found for these two-loop EHL's C^ in (F) (23J [HI ESJ [26j E3 EH ES] and 

42a (-F) [301 El EH EH [29], all at the two-parameter integral level. So 
far these integrals have withstood attempts at analytical evaluation. How- 
ever, they have been used to compute the two-loop weak field expansion 
coefficients to high orders |26|, [27] [22] , and to derive closed- form expressions 
for these coefficients in the purely electric/magnetic case [31] . 

The Schwinger representations (|1.5|) . (|1.6|) also generalize to the two-loop 
level [32[ I24j. although now with a certain prefactor function Kk(eE/m 2 ) in 
front of the A;th Schwinger exponential; those functions are known explicitly 
only at leading orders in the weak-field limit. 

However, it turns out that the electric or magnetic backgrounds are not 
the simplest ones in this context. Computationally, the most favorable case 
is the one of a (euclidean) self-dual ('SD') field, defined by F = F, which 
has the consequence that 

F 2 = -f 2 t (1.10) 

For a real value of the parameter /, the SD effective Lagrangian has proper- 
ties similar to the magnetic EHL, for imaginary / similar to the electric one. 
In this SD case, even at two loops it is possible to obtain explicit formulas 
for the EHL's [H1CE2]: 
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Here we have defined the convenient dimensionless parameter 
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m 

as well as the function 
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where ip is the digamma function ip(x) = T'(x)/T(x). Thus in the self- 
dual case the study of the weak-field expansions and the construction of the 
imaginary parts involve only some well-known properties of the digamma 
function, making it possible to do everything much more completely and ex- 
plicitly than for the electric or magnetic case [12] . In particular, for this case 
it has been verified that certain complications related to the non-uniqueness 
of /, which could invalidate the relation between (II. 8p and (II. 9j) . such as the 
appearance of additional poles or cuts in the complex g plane, do not occur 
in QED, or at least not yet at the two-loop level. 

Moreover, although a self-dual field cannot be realized in Minkowski 
space, the corresponding EHL still contains information on the N - photon 
amplitudes, namely on their "all +" component in the helicity decomposition 

is]. 

Beyond the two-loop level, to the best of our knowledge the only result 
on EHL's in the literature is the following formula, proposed in 1982 by 
Affleck, Alvarez, and Manton |34j as an all-loop generalization of the leading 
Schwinger exponential for scalar QED (jl.5p : 
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This formula (called 'AAM formula' in the following) is highly remarkable 
for various reasons. Despite of its simplicity it is a true all-loop result; 
the rhs receives contributions from an infinite set of Feynman diagrams of 
arbitrary loop order, as sketched in fig. [TJ 
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Figure 1: Diagrams contributing to ImC 



(all— loop) 
seal 



(E) in the weak-field limit. 



Moreover, the mass appearing in (jl . 15[) is argued to be still the physical 
renormalized mass, which means that the above figure should strictly speak- 
ing include also the mass renormalization counter diagrams which appear in 
EHL calculations starting from two loops. 

The derivation given in [3l] is very simple, if formal. Based on a station- 
ary path approximation of Feynman's worldline path integral representation 
[35] of 

>Cscai (-£■); h actually uses only a one- loop semiclassical trajectory, and 
arguments that this trajectory remains valid in the presence of virtual pho- 
ton insertions. This also implies that non-quenched diagrams do not con- 
tribute in the limit (|1.15p . which is why we have shown only the quenched 
ones in fig. [TJ 

Although the derivation of (|1.15p in [M] cannot be considered rigorous, 
an independent heuristic derivation of (|1.15p . as well as extension to the 
spinor QED case (with the same factor of e a7r ) was given by Lebedev and 
Ritus [32] through the consideration of higher-order corrections to the pair 
creation energy in the vacuum tunneling picture. At the two-loop level, 
(|1.15p and its spinor QED extension state that 
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which has been verified by a direct calculation of the EHL [22J (for the spinor 
QED case), and also been extended to the case of a self-dual field 112]. 

Now, writing the AAM formula (|1.15p as (hereafter, unless indicated 
otherwise, we refer to the scalar and spinor QED cases simultaneously) 

OO 

I m £(ali-ioo P )^ = J2~kn£ (l) (E) E ~° Im£W(£) e a7T (1.17) 
1=1 

it states that an all-loop summation has produced the convergent factor e Q7r , 
clearly an observation similar in vein to Cvitanovic's. Moreover, at a formal 
level it is not difficult to transfer this loop summation factor from Im£(E) 
to the QED photon amplitudes \12\ [14] . Consider the weak field expansion 
of the Z-loop contribution to the electric EHL: 

c®w = fyv)gf (i.i8) 

n=2 

Using the Borel dispersion relations discussed above, f 1 1 . 1 5 1) can be shown 
[22] [T2] to imply that, at any loop order I, the weak field expansion coeffi- 
cients have the same leading asymptotic growth as we found above for the 
one-loop case, that is 

S(n) n ~°° c« 7T- 2n T(2n - 2) (1.19) 

where the constant relates directly to the prefactor of the corresponding 
leading Schwinger exponential in the weak field limit: 

ImC^(E) c«e-^ (1.20) 

We would now like to convert (jl.20p into an equation for photon amplitudes, 
using the above-mentioned correspondence between the weak field expansion 
coefficients and the EH limit of the photon amplitudes. Here we encounter 
the problem that this relation cannot be applied to a purely electric field; we 
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will therefore now switch to the SD case, and, on the strength of the two-loop 
results of I12j. assume that (|1.17p holds for the SD case unchanged. 

Now, the SD EHL relates to the 'all +' N - photon amplitudes, and for 
those the whole kinematic structure in the EH limit can, independently of 
the loop order, be absorbed into the invariant x% defined in (jl.3p [13J. One 
can thus eliminate this kinematic factor by dividing the I - loop amplitude 
by the one-loop one. Expanding (|1.15j) in a and combining it with (|1.20|) 
and (|1.19p one then arrives at a formula for this ratio of amplitudes in the 
limit of large photon number [T2] , 



r r(0[fc 1 , £ +;...;fc JV , g +] c O(n) _ (an) 1 ' 1 

hmN ^m kl ,ei;...;k N ,e%] ~ hm ™° ~ J^l 

(1.21) 

If we could now sum both sides over I and interchange the sum and limit, we 
could reconstruct the e a7T factor, and conclude that the perturbation series 
for the iV - photon amplitudes, if only for sufficiently large N, in the low 
energy limit and for this polarization component, has infinite convergence 
radius. But this is too good to be true, since so far we have nowhere made a 
distinction between quenched and unquenched contributions to the photon 
amplitudes, and convergence of the whole perturbation series can certainly 
be excluded. 

However, as was noted in |14j this distinction comes in naturally if one 
takes into account that, as mentioned above, in the path integral derivation 
of (|1.15p in [34] the rhs comes entirely from the quenched sector; all non- 
quenched contributions are suppressed in the weak field limit. And since 
(switching back to the usual Feynman diagram picture) the importance of 
non-quenched diagrams is growing with increasing loop order, it is natural 
to assume that their inclusion will slow down the convergence towards the 
asymptotic limit with increasing I, sufficiently to invalidate the above naive 
interchange of limits. On the other hand, there is no obvious reason to expect 
such a slowing down of convergence inside the quenched sector, which led 
to the prediction [14j that Cvitanovic's "quenched convergence" will indeed 
be found to hold true for the photon amplitudes. 

To further corroborate this prediction, one should now calculate the 
EHL's and their weak field expansion coefficients c^ 3 \n) at the three-loop 
level. This would already allow one to perform some quite nontrivial checks 
on the above chain of reasoning, namely: 
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1. The asymptotic relation 
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and its independence of spin. 

2. The absence of a slow-down of the convergence of the ratio (jl.22p as 
compared to the corresponding two-loop to one-loop ratio. 

3. The asymptotic subdominance of the non-quenched part of the EHL, 
which first appears at three loops. 

However, a calculation of any three-loop EHL, be it in Scalar or Spinor 
QED, for an electric or self-dual field, presently still poses an enormous 
computational challenge. 

Now, it is well-known that the structure of the EHL's and associated 
Schwinger exponentials at one-loop is essentially independent of the space- 
time dimension [36j [371 Ell EH]. In particular, the Schwinger formulas for 
1+1 dimensional QED are different from (|1.5j) . (|1.6p only by normalization 
factors: 
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More surprisingly, in 2006 Dunne and Krasnansky [401141] found that, for the 
scalar QED case, even at the two-loop level there persists a strong similarity 
between the 2D EHL and the 4D self-dual EHL. For the former, they find 



= -£| (1.25) 

2 

where a := 2-^-r is our definition of the fine structure constant in two 
dimensions, and 
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:= - (^(« + j) - ln ( K )) = " 2 ^( 2 «) + ln(4«) 

(1.26) 

Here the parameter / is defined by (see app. [A]for our 2D QED conventions) 



F - ( -/ o ) < L27 > 

(note that this definition is consistent with (|1.1U|) ). The formal similarity 
between ([L25]) and the 4D SD EHL (fTTTTT) has led us to consider 2D QED as 
a toy model for studying the above asymptotic predictions. This will make 
sense only, of course, if the AAM formula fjl.28j) can be generalized to the 
2D case. As we will show below, this generalization is 



2 ~ r„2„2 
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As in the 4D case, we shall assume that spin does not play a role in this limit, 
so that (jl.28p holds for spinor QED unchanged (let us also mention here that 
Loskutov et al. [42] have found an all-loop exponentiation formula similar 
to (jl.28p for the strong-field asymptotics of the electron mass operator in a 
two-dimensional approximation of four-dimensional QED). 

In the same way as in the 4D case, one can then use Borel analysis to 
derive from (|1.28j) a formula for the limits of ratios of I - loop to one - loop 
coefficients: 



ciUn) (an 2 ) 1 - 1 . . 

lim^oo 2DK ' = \ [ 1.29 
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Here our definition of the expansion coefficients in 2D is 



n=l 

Note that, due to the factor of k 2 appearing in the exponent on the rhs 
of (|1.28|) . the 2D formula (|1.29|) involves also a shift in the argument of 
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c 2Z)( n )' ^ e leading asymptotic growth of the coefficients increases with the 
loop order. This actually simplifies matters, since it implies that, for our 
asymptotic purposes, there is no need to perform mass renormalization: 
Mass renormalization contributions to an EHL at the I - loop level are of 
the form 

5m {h) 5m (h) ■ ■ ■ 5m^—- C {V) (1.31) 

with some V < I and I — V = Xa=i h> where m^ li ' denotes an ^ - loop mass 
counterterm. As we will show below, Co^(n) grows like T(2n — 1), and (jl.29p 

then predicts that c^ D (n) grows like T{2n + 21 — 3). It is then easily seen 
that the expansion coefficients of a term of the structure (jl.3ip at I loops 
can grow at most like ~ T(2n + 21 — 4), and thus are subdominant. 

This is very different from the 4D case. Here the leading asymptotic 
growth of the expansion coefficients of the unrenormalized electric EHL's 

42o2n is like r ( 2n + 1 ~ 3 ) ( F ( 2n + 1 ~ 2 ) for the SD case )' that is ' the 
argument of the T - function jumps only by one per loop order |14| . Thus 

mass renormalization terms do contribute at the leading asymptotic level, 
and are even crucial to make the AAM formula work: At two loops, it has 
been shown [22 1 [T2] that, precisely when the mass is taken to be the phys- 
ical renormalized mass, the leading asymptotic terms cancel between the 
main contribution and the one from mass renormalization, which reduces 
the asymptotic growth at two loops to make it the same as at one-loop, 
as is implicit in the AAM formula f 1 1 . 1 5 1) . And for this formula to hold, 
these cancellations between main and counterterms must not only persist, 
but become increasingly extensive at higher loop orders. Thus the AAM 
formula already predicts that the asymptotic behaviour of QED must de- 
pend crucially on whether mass renormalization is done physically or just 
generically. 

In section [2] we will first retrace the original derivation of the AAM 
formula and its generalization to spinor QED |43j . and then present our 
generalization to the 2D QED case. In sections [3] and [4] we compute the one- 
and two- loop EHL's in 2D spinor QED. In section [5] we use these results, 
together with the formulas previously obtained in [401 141] for scalar QED, 
to verify that the asymptotic predictions derived from our generalized AAM 
formula hold at two loop order in both theories. Section [6] gives a summary. 
The second part of this series |44| will be devoted to the calculation of the 
three-loop EHL's in scalar and spinor 2D QED, and to the verification of 
the 2D analogue of the three-loop predictions discussed above. 
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Some preliminary results of this work have been presented in [45]. 



2 Predictions of the AAM formalism 



In this chapter, we will first retrace the derivation of the AAM formula 
(|1.15p . and then use their method to obtain also the generalization of this 
formula to the 2D case, eq. (|1.28p . 

The seminal work of those authors concerned scalar QED in four dimen- 
sions, therefore let us start with this case, beginning at the one-loop level. 
The (euclidean) one-loop effective action for scalar QED can be written in 
the following way |35| : 



seal 
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Here m is the mass of the scalar particle, and the functional integral J T>x is 
over all closed spacetime paths x^(t) which are periodic in the proper-time 
parameter r, with period T. Rescaling t = Tu, the effective action may be 
expressed as 



r sca i[A] = 
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After this rescaling we can perform the proper-time integral using the method 
of steepest descent. The T integral has a stationary point at 
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mi 



dux 2 ^ \ 



(2.3) 



(2.4) 



(2.5) 



12 



The functional integral remaining in the effective action expression ([27 
may be approximated by a further, functional, stationary phase approxima- 
tion. The new worldline "action" , 




-1 r-l 

S e s = m \l I dux 2 + ie duA-x (2-6) 
is stationary if the path x a {u) satisfies 

m Xfl = = ieF^Xy (2.7) 
yj Jldux 2 

A periodic solution x^(u) to (|2.7p is called a "worldline instanton". Con- 
tracting (|2.7|) with Xu shows that for such an instanton 

x 2 = constant = a 2 (2-8) 

Generally, the existence of a worldline instanton for a background A leads to 
an imaginary part in the effective action r sca j[A], and the leading behavior 
is 

Imr scal L4] ~e~ s ° (2.9) 

where So is the worldline action (|2.6p evaluated on the worldline instanton. 

For a constant electric background of magnitude E, pointing in the z 
direction, the Euclidean gauge field is A^^x^) = —iEx^. The instanton 
equation (|2,7p for this case can be easily solved, and the solutions are simply 
circles in the z — t plane of radius 



x\{u) = — cos(2A;7rn) , x\{u) = sin(2£;7ru) (2.10) 

(with x\ 2 kept constant). The integer k G Z + counts the number of times 
the closed path is traversed, and the instanton action (|2.6j) becomes 



S :=S eS [x»] = 2k — -k— = k— (2.11) 



Thus in the large mass approximation (|2.5p the contribution of the instan- 
ton with winding number k reproduces the exponent of the kth. term of 
Schwinger's formula (|1.5|) . 
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Proceeding to the multiloop case, Feynman's formula (|2.ip can be easily 
modified to take into account the effect of multiple exchanges of photons in 
the scalar loop. This requires only the addition of the following interaction 
term Si to the worldline action, 



e2 f T dTl [ T dT2 _*(n)- X (r 2) 



8^J >(ti)-*(t 2 ))2 

(2.12) 

After this addition in (|2.ip . r sca i[A] turns into the "quenched" effective ac- 
tion, corresponding to the set of all Feynman diagrams with a single scalar 
loop (but any number of interactions with the external field and internal 
photon exchanges). Now, Affleck et al. [34J argue that the instanton so- 
lutions (|2.10p remain stationary even in the presence of this addition, so 
that its effect in the large mass limit is only a modification of the stationary 
action. Evaluation of Si on the leading k = 1 instanton yields simply 



Si[x{u)} = -air (2.13) 

Combining (|2.9p . (|2.1ip . (|2.13p one obtains the large mass approximation for 
the quenched multiloop scalar EHL (|1.15p . 
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Affleck et al. then proceed to show that the contributions to the full effective 
action involving more than one scalar loop are suppressed with respect to 
the quenched one in this large mass limit, so that f|l . 15[) actually holds even 
for the full effective action. Even more remarkably, they argue that the 
asymptotic formula (11 . 151) already takes into account all the effects of mass 
renormalization which come into play for this effective action starting with 
the two-loop level. 

The approach of Affleck et al. was generalized in [13] to the case of spinor 
QED. The path integral representation of the one- loop effective action due 
to a spin half particle in the loop differs from (|2.ip only by a global factor 
of — \ and the insertion of the following 'spin factor' 5[x,j4] under the path 
integral [16] . 

S[x,A]=ti r V e^/o^^M-)) ( 2 .14) 
14 



Here trp denotes the Dirac trace and V the path ordering operator. The 
presence of this spin factor term in the stationary path approximation is 
again taken into account simply by evaluating this term on the stationary 
path (|2.10p . The result turns out to be a simple global factor [4~3] , 



4cos(™) =4(-l) n (2.15) 

Taking the global factor of — \ into account this correctly reproduces the 
difference between the one-loop Schwinger formulas for scalar and spinor 
QED (|1.5p . (|1.6p . Proceeding to spinor QED at the multiloop level, here there 
seems to be no simple way to directly calculate the effect of the spin term 
in the stationary path approximation. However, as was discussed already in 
the introduction the physical arguments of [32] as well as explicit two-loop 
calculations strongly suggest that the AAM formula (|1.15p holds also for the 
spinor QED case (up to a global factor of 1/2). 

We will now generalize the AAM formula to the case of QED in two 
dimensions. First, we note that the path integral representation (|2.ip for the 
scalar QED effective action is independent of dimension; the only explicit 
appearance of the dimension D is in the normalization of the free path 
integral, which is 



x(T)=x(0) 



T>x exp 



(4ttT)' 



-D/2 



(2.16) 



(with the conventions of [29] ) . This does not affect any of the manipulations 
which we performed at the one-loop level, including the form of the worldline 
instanton (I2.10p . so that the asymptotic estimate (I2.9P remains valid in 
D = 2 with the same formula (|2.1ip for So. 

Things are different at the multiloop level. Here it must be remembered 
that the worldline insertion term (|2.12p actually involves the photon prop- 
agator in D = 4 and in Feynman gauge. For general D and in a general 
covariant gauge this term reads 



Si{D) 



e 2 1 

2 ■ D 



4-7T 2 J 



T 



dr a / dr b < 



(x a - x b y 



i _ a )r( D ) Xa ' ( Xa - Xb ^ Xa 



V 2 



x b ) • x b 



(2.17) 



(x a - x b y 
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where a = 1 corresponds to Feynman gauge. Since the first term in braces 
in (|2.17p becomes singular in D = 2, in this case instead of Feynman gauge 
it is more convenient to choose the gauge a = — 1, leading to 



5 f=2 = £i r dTa f T dn ^i^i x b)( x « x b)-i b (218) 

47r Jo Jo (sc« - x b y 

The evaluation of on the constant field worldline instanton (|2.10p then 
yields 



S?=Mu)] = ^ = 2vr(m) 2 ^ = -avrV (2.19) 
This brings us to our 2D generalization of the AAM formula, eq. (|1.28D . 



3 1 loop Euler-Heisenberg Lagrangian 



The (renormalized) one-loop Euler-Heisenberg Lagrangians in 2D scalar and 
spinor QED are, in the standard proper-time representation (see, e.g., [36]), 



et/) - -gf fe-(^)-i) (3.2) 

It is convenient to observe that the scalar EHL can be written in terms of 
the spinor EHL as [17] 



spin 



(/) 



2/£l,(//2) 



spin v 



(3-3) 



The integrals (|3.ip . (|3.2p have various equivalent closed-form representations; 
the most suitable one for our purposes is 



m 1 

47T ft 



lnr(K) — «(ln k — 1) 



2 V 2vr ; 



(3.4) 
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4 2 loop Euler-Heisenberg Lagrangian 



We now calculate the two-loop Euler-Heisenberg Lagrangian. In momentum 
space, it is given by 



After plugging in ()A.8|) , we note that due to (jA.lOp the terms involving a ■ p 
do not contribute under the Dirac trace, which thus becomes simply 



tr [e^V^^ 2 '^] = tr [ e^ 2 e^V^] = 4cosh(z-z') (4.2) 
This leaves us with 

v Jq Jq cosnzcosnz' 

d 2 p f d 2 p' 1 e --L,(tanhzp 2 +tanhz'p' 2 ) 



(2n) 2 J (2tt) 2 (p - p 1 ) 2 
We exponentiate the photon propagator, 



1 



oo 

,'12 



dXe-^P'y (4.4) 



(p ~ P') 2 Jo 
The momentum integrals can now be done, yielding 



(fp I d y e -^(P-P') 2 -^(tanhzp 2 +tanhz'p' 2 ) _ ^(e/) 2 



tanhztanhz' + A ef(tanhz + tanhz') 

(4.5) 



The A - integral has a logarithmic infrared divergence which we regulate 
with a cutoff Ao, / °° dX — > dX. Using 

' ' (4-6) 



cosh z cosh z'(tanhz + tanhz') sinh(z + z' 

the result becomes, in the limit Age/ ^> 1, 
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^spin(/) 



IM_ r dT r dT , c -mi(T + T ^™Hz-z') 

8vr 2 Jo J smh(z + z') 



ln(e/A ) + lnf sinh ( z + z ') 
V sinh z sinh z' 

(4 



For the calculation of the integrals, it is convenient to change to the "world- 
line" variables 1291 . 



T 



uT 

{l-u)f 



oo roo 



dT \ dT' 

n Jo 



oo rl 

dff / du 
o Jo 



(4.8) 



There are two different u integrals, both elementary: 



lilt cosh [(1 — 2u)Z~\ 



cfot cosh [(1 — 2u)Z] lnsinh(uZ) 



sinh(Z) 



- — cosh Z + — fin sinh Z ) 

9 v. V o/ 



2' 
(4.9) 



where Z = e/T. After some rearrangements, the result can be written as 



£ {2) . (f) 



m 2 e 2 1 rQO 



Svr 2 e/y 



dZe" 



-2kZ 



Z(cothZ — — ) — In sinh Z + In Z 



+ 



8?T 2 L 



ln(A m 2 ) +2 + 7 



(4.10) 



(7 is the Euler-Mascheroni constant). The remaining integrals can be easily 
reduced to the standard integral [48J 



f°° 1 
/ dZ e~ 2KZ (cothZ - - 
Jo Z 



(4.11) 



Putting the pieces together, our final result for the (unrenormalized) 2-loop 
spinor EHL in 2D is 



18 



42» = ^^[-e , (^) + ln(A m 2 ) + 7 + 2] (4.12) 

Thus, contrary to the scalar QED case, the two-loop 2D spinor QED 
EHL turns out to be significantly simpler than the corresponding self - dual 
4D one f)l . 12|) . In fact, it has the curious property that, up to the vacuum 
part, one can even write it in terms of derivatives of the one- loop EHL (13.4[) : 



43„w = -!('" 2 i) 2£ S>> (4.i3) 

Note that, as explained in the introduction, there is no need for us to perform 
the mass renormalization of C^ in . 

5 Weak field expansions 

Finally, we will now work out the weak field expansions of the one and 
two-loop 2D EHL's, for both the scalar and spinor cases, and verify the 
asymptotic prediction (jl.29p derived from the AAM formalism. 

Starting with the spinor case, for working out the one-loop EHL (|3.4h 
we need the large x expansion of lnr(a;), which is 



1 °° p 

lnr(.) ~ .(In. - 1) -iln(J) + £ ^TTf^ (5.1) 

n=l v ' 

This yields (the subscript '2D' will be omitted in this chapter) 



<&<"> = ^^T) < 5 ' 2 > 

(n > 1). At the two- loop level, we get from (15. ip that 



~ E^r^" 2 " (5.3) 

n=l 

Thus (|4.12p gives (omitting the cutoff-dependent vacuum energy term) 
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<&(») = i-V n+1 l 2 -^B 2n (5.4) 
(n > 1). Using Euler's formula 

B 2n = (-l) n+1 2(2vr)- 2n (2n)!C(2n) (5.5) 
and lirn n ,_ >00 £(n) = 1 we find 



(2) / s «r(2n+ 1) 



spinv > 4 ( 27r )2 



71 



(5.6) 



so that 



,. SspinVV ~ 2 / r >-7\ 

lim -7— = air (5.7J 



<Spin(« + 1 ) 

in agreement with (jl.29|) . 

Proceeding to the scalar QED case, at one-loop we get, from (|3.3p and 

63, 

(n > 1). Comparing with (|5.2|) . we see that c£^(n) ~ c lp| n ( n ) f° r large 
n. At two loops, we need the expansion of the function £,2d(k), defined in 
(fT2"5]) , which is 



DO 



n=l 

where we have defined B n := (2 1 "™ - l)B n . The expansion of (jl.25p then 
gives 
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n-1 ^ 



B2m B2n-2m 



i 2 i(») - £(-ir(E 



4m 4(n — m) 




(5.10) 



m=l 



(n > 2). The structure of these coefficients is thus very similar to the one of 
the expansion coefficients of the two-loop EHL's for a self-dual field in 4D 
[T2] . In [12] also a method was developed to compute the asymptotic 
expansion of folded sums of Bernoulli numbers of the type appearing in the 
first term in brackets in f|5. lOj) ; applying the same technique to the case at 
hand one can show that this term is subdominant in the large n limit. The 
leading order contribution comes from the second term, and yields 



Comparing with (|5.4|) . we see that the two-loop coefficients c scal (n) and 

c spin( n ) a ^ so a g ree f° r large n. This completes our check of (|5.7p for the 
scalar case. 

6 Conclusions 

To summarize, in the first part of this series we have, generalizing the work 
of Affleck et al., obtained an all- loop formula for the imaginary part of the 
Euler-Heisenberg lagrangian in 1+1 QED, in the limit of a weak field but 
at arbitrary coupling. We have performed the first calculation of the Euler- 
Heisenberg lagrangian in 2D spinor QED at two loops, and verified that the 
asymptotic behaviour of its weak field expansion coefficients agrees with the 
prediction of this generalized AAM formula. We have also performed the 
same check for the scalar QED case, using the two-loop result of Dunne and 
Krasnansky. Our findings clearly demonstrate that 2D QED is sufficiently 
close to the 4D case to suggest that it may hold generic information on 
the asymptotic behaviour of amplitudes in QED in general, but at the same 
time holds more promise for the explicit study of multiloop Euler-Heisenberg 
lagrangians. A particularly interesting common aspect of our two-loop 2D 
QED results and the two- loop SD EHL calculations of [T2] is, that, in 
all cases where an explicit formula has been obtained for a two-loop EHL, 
the leading asymptotic growth of its expansion coefficients has turned out 




aT(2n + 1) 
4 (2vr) 2 ™ 



(5.11) 
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to involve only the digamma function, and only linearly. We take this as 
further evidence that the study of the iV-photon amplitudes through the 
QED effective Lagrangian may ultimately provide a window to high orders 
in perturbation theory. 
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A Conventions and formulas for Euclidean 1+1 
QED 



Dirac equation: 








(A.l) 



Free electron propagator: 



1 



— ip 1 + m 



(A.2) 



ip 1 + m 



p2 _|_ m 2 



Photon propagator in Feynman gauge: 



5, 




(A.3) 



Vertex: 



(A.4) 



Constant field: 



F 




(A.5) 



Fock-Schwinger gauge: 
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Electron propagator in a constant field in Fock-Schwinger gauge: 



(A.6) 



[(T IJt (df i -ieA ft ) + m]g(x-x') = 5(x - x') (A. 7) 



9(P) = r^e- T(m2+ ^ p2) ^T-f-e CT32 -^) (A-8) 
Jo cosh z V cosh z J 

9{x) = Iwf eipX9{p) (A - 9) 

1 f°° dT m 2 T z z J ( a 3 z 1 z a 

= — / — e" m 1 — - — e'wSff me ff ^ 

47r 7 T sinh z V 2T sinh z J 

(z = efT). 

One of the motivations for considering the Euler-Heisenberg Lagrangian 
in 1 + 1 dimensions is that substantive simplifications can be expected for 
higher loop calculations in Feynman gauge due to the fact that 

er M 0Vi • • •0V 2 „+i0> = (A. 10) 
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